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We study the interplay between large-spin, spin-orbit coupling, and superfluidity for bosons in a 
two dimensional optical lattice, focusing on the spin-1 spin-orbit coupled system recently realized 
at the Joint Quantum Institute [Campbell et. ah, arXiv: 1501.05984 . We find a rich quantum phase 
diagram, where, in addition to the conventional phases —superfluid and insulator— contained in the 
spin-1 Bose-Hubbard model, there are new lattice symmetry breaking phases. For weak interactions, 
the interplay between two length scales, the lattice momentum and the spin-orbit wave-vector induce 
a phase transition from a uniform superfiuid to a phase where bosons simultaneously condense at the 
center and edge of the Brillouin zone at a non-zero spin-orbit strength. This state is characterized 
by spin density wave order, which arises from the spin-1 nature of the system. Interactions suppress 
spin density wave order, and favor a superfiuid only at the Brillouin zone edge. This state has 
spatially oscillating mean field order parameters, but a homogeneous density. We show that the 
spin density wave superfiuid phase survives in a two dimensional harmonic trap, and thus establish 
that our results are directly applicable to experiments on ®^Rb, ^Li, and 

PACS numbers: 67.85.Bc, 67.85.Jk, 67.85.Fg, 67.85.Hj 


The interplay between spin-orbit coupling, lattice 
physics, and interactions is fundamental to many areas 
of condensed matter and materials physics from topo¬ 
logical insulators^, lattice analogues of quantum Hall ef¬ 
fects^, spin liquids^, and spintronics based devices^. Re¬ 
cent progress with ultra-cold atoms and molecules has al¬ 
lowed us to “engineer” spin-orbit coupling and competing 
interactions in systems that have no solid state counter¬ 
parts. This enables the study of fundamental questions 
such as the fate of superfluidity/superconductivity in the 
presence of single particle degeneracies^-—, the physics of 
bosons and fermions with large spini^ii^, quantum Hall 
effects in bosonic systemsi^, and the physics of long range 
inter act ionsi^Ti^S. In many cases, the very paradigms to 
qualitatively think about these strongly correlated sys¬ 
tems is only now being developed. 

We theoretically study the interplay of spin-orbit cou¬ 
pling, large spin, and interactions by calculating the 
ground state phase diagram of a spin-1 spin-orbit cou¬ 
pled Bose gas in an optical lattice. In the continuum, 
spin-orbit coupling (SOC) gives rise to a single particle 
energy dispersion with multiple minima, each with a dif¬ 
ferent spin character. For weak interactions, and large 
on-site occupation number, bosons condense into mul¬ 
tiple minima simultaneously, producing a ground state 
with a spatially oscillatory order parameter and density 
modulations21“i2i, i.e. a stripe superfiuid (SF). In the 
spin-1 system, the spin also develops a spatial texture 
through a spin density wave (SDW) SF— . This SDW 
SF has yet to be observed experimentally in any cold 
atom system. 

As we show here, the physics in a lattice is strikingly 
different. In addition to inducing strong correlations by 
suppressing inter-particle tunneling, the lattice also sup¬ 
presses any SF formation at incommensurate momentum. 



FIG. 1: (a) Schematic representation of the model on a square 
optical lattice. Hopping in the x-direction is spin depen¬ 
dent carrying a SOC driven phase, while hopping in the y- 
direction is spin-independent with no phase, (b) Mean field 
phase diagram of the spin-1 spin orbit coupled Bose-Hubbard 
model with U 2 — —0.25t/o, TLr — O.lHo and 0 — 'k. We 
find three distinct spin-orbit driven phases: a S-SF, a Q-SF, 
and density (A) tuned Mott phases at strong coupling (the 
phases are described in the text). Interactions suppress the 
SDW order leading to a second order transition from the S- 
SF into the Q-SF phase, where the SF is now ferromagnetic 
and condensed at q = (tt, 0). Complete translation symmetry 
is restored at a second transition to the ferromagnetic Mott 
phase. Our conclusions remain qualitatively valid even for 
weak U 2 — —0.005f/o, which is the case for ®^Rb. 


As a result, instead of a stripe superfiuid modulating at 
the spin-orbit wave-vector for arbitrarily weak spin-orbit 
coupling, the lattice induces a second order phase transi¬ 
tion between a homogeneous ferromagnetic (FM) SF to a 
stripe SDW superfiuid (S-SF) at a strong, non-zero SOC 
strength. Remarkably, this stripe superfiuid involves con¬ 
densation at the center and the edge the Brillouin zone, 
and not at momenta (necessarily) commensurate with the 
spin-orbit wave vector. This phase transition is a direct 
result of the interplay between two length scales, namely 
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the SOC wave-vector where the bosons prefer to condense 
in momentum space, and the lattice momentum, which 
describes where the bosons are pinned in real space. 

Approaching the Mott limit simultaneously suppresses 
the superfluid and stripe order— and favors conden¬ 
sation at a single momentum, due to the competition of 
the various energy scales in the problem. A key question 
then becomes whether SDW order and superfluidity are 
simultaneously destroyed at the Mott transition, or are 
these are two distinct quantum phase transitions? We 
find that the SDW superfluid is destroyed well before the 
Mott boundary at a continuous transition, giving rise to 
an interaction driven SF phase where condensation oc¬ 
curs at q = (tt, 0 ) on the Brillioun zone-edge. This novel 
superfluid has real space oscillations in the superfluid or¬ 
der parameter, but not in the density. Our predictions 
can be readily tested in experiments by generalizing the 
setup of Campbell et al^ to a two dimensional optical 
lattice. 

Spin -1 Bose-Hubbard model with SOC — We consider 
the spin-1 Bose-Hubbard model^^^ in two dimensions, 
in the presence of a one-dimensional spin-orbit coupling 
H = W + Y.iHl 

+ ( 1 ) 

i'i'J) 

Hi = ^ni{ni - 1 ) + ^ (Sf - 2ni) + ^IrSI + ViUi, 

where describes the creation of a spin-1 boson in the 
ruz = a state, repeated Greek indices are summed over, 
and denotes a sum over nearest neighbors. We in¬ 
troduce the density Ui = a\^aio,^ the spin = aJ^F^^ai/3 
(where F is a vector of the spin-1 matrices), the spin-orbit 
coupling [via the spin dependent hopping see Fig. [T] (a)] 
and = 1, the 3 x 3 identity matrix, 

the density-density interaction I/q, and the spin-spin in¬ 
teraction 1/2- We study repulsive interactions Uq > 0 , 
and a ferromagnetic spin-spin interaction I/2 < 0, as for 
^^Rb, ^Li, and (Ref. 1 ^ . We work at a fixed chem¬ 
ical potential /i, which enters the calculations through 
Vi = Wi — where Wi is the harmonic trapping poten¬ 
tial. We include a Zeeman coupling (i.e. Raman coupling 
strength) CIr perpendicular to Sx, which is essential to 
ensure the SOC cannot be simply gauged away by a lo¬ 
cal unitary rotation. Importantly, the physics we discuss 
does not require strong spin-orbit coupling, and thus will 
not suffer from heating losses. 

We use an inhomogenous Gutzwiller mean 
field theory (GMFT)^, with a wavefunc- 
tion ansatz 1^) = Where \i>i) = 

is the Wave 

function at site i, and we determine the 
variationally, in an unbiased fashion allowing them to 
be different at each lattice site, while truncating each 
local Hamiltonian to a bosonic basis In the absence 
of SOC, the GMFT qualitatively describes polar and 
ferromagnetioi^ii^ superfluid phases correct l}^”—, and 



FIG. 2: Evolution of the ferromagnetic superfluid as a func¬ 
tion of the spin-orbit angle 0 for Lx — Ly — 20 for LLr — 
O.R/o, and t — 0.21/o [for panels (a), (b) and (d) U 2 — 
—0. 251 / 0 ]. Absolute value of the real and imaginary (inset) 
parts of the superfluid order parameters for mz — 1 (squares), 
0 (circles), and —1 (triangles) at a momentum q = (0, 0) (red) 
and q = (tt, 0) (blue) as a function of 0 . (b) The momentum 
dependence of |(a(q))| = y^|ai(q)|2 + |ao(q)|^ + |a_i(q)|2 for 
= 0, displaying that the weight of the mode at q = (0,0) 
gets transferred to q = (7r,0) at the critical angle Oc- (c) 
The spin density wave (SDW) order parameter as a function 
of 0 showing a continuous transition into the S-SF phase for 
I/ 2 /I /0 = —0.005 and —0.250. (d) The spin {S^{x,y)) (in 

units of the average density) in real space, displaying a trans- 
lationally invariant solution in the ^-direction and a distinct 
SDW pattern. 


also captures the presence of each Mott lobe. Here, we 
explore how this phase diagram is modified in the pres¬ 
ence of the spin-orbit coupling, and in particular, how 
the superfluid and spin density wave order parameters 
evolve close to the Mott transition. It is important to 
remark that quantum Monte Carlo^, which is often 
exact for bosonic systems incurs a sign problem from the 
complex hopping. Variational mean-field theories of the 
form we consider are thus essential to making progress. 

We first calculate with no trap (Wi = 0 ) and 

periodic boundary conditions (Figs [U [2l and [ 3 ]). The 

explicit form of (V’l {'4>\aia\^), and (^/'|o|„|?/)) 

can be found in Ref. l 38 l. while allowing the variational 
parameters to be site dependent. The hopping Hamilto- 
nian becomes +-Ria(«*«), 

+ (7^y)*l (a,v3), and Ria = 




where Ria = 


where ^ 


{3)i 


denotes a 


sum over the nearest neighbors of site i. We diagonalize 
each site of the resulting local Hamiltonian, which feels 
the neighboring condensate order parameter (through 
Rio) iteratively until the lowest energy state is reached. 
To avoid getting trapped in local energy minima, we solve 
the mean field equations starting from different random 
initial conditions (for the (a^c^)), and find the solution 
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with the lowest energy. We consider small Qr = O.lUo^ 
to avoid simply polarizing the spinor gas^. In this work 
we focus primarily on the deep lattice limit (t ^ Uq) 
and low particle number (considering bosonic Mott lobes 
up to A/" = 3 ), and therefore can restrict the calcula¬ 
tions to a small number of bosons per site. As we show 
below, solving the fully two dimensional inhomogenous 
problem yields translationally invariant solutions in the 
^-direction, and a two site unit cell in the x-direction 
[see Fig. [ 2 ][d)]. Thus in constructing the phase diagram 
in Fig. (TJ we consider translationally invariant solutions 
in the ^-direction. 

Striped Superfluid and Spin Density Wave Phase — At 
weak coupling, the physics is dominated by the interplay 
of lattice effects and SOC, which we explore in Fig. [2j 
We fix U2 = — 0 .251/0, / = O.2I/0 for lattice sizes = 
Ly = 20 (setting the lattice spacing to unity), and choose 
a small value of the chemical potential /i = 0 . 5 //o, such 
that a = 3 basis is sufficient on each site. We remark 
that this value of I/2 is experimentally realistic as ^Li has 
U2/U0 ^ — 0.5 (Ref. [H). We take the Fourier transform 
of the bosonic operators to determine their structure in 
momentum space via adq) = L~^Ly^ For 

6> = 0, the ground state is a ferromagnetic superfluid 
with (a_i(q = 0)) 7^ 0 (where the term containing CIr 
in Eq. ([ 2 ]) selected the ordering direction), as shown in 
Fig. [ 2 ][a) and A spiral magnetic state 

develops with increasing 0 [Fig. [ 2 ] (a)] as (ai(q = 0 )) 7^ 0 
for small 0 . 

As 0 increases further, we find a continuous quantum 
phase transition at 0 c{t^U 2 = — 0 . 25 //o, f^i?) ~ 2.2 into 
a SDW superfluid phase (S-SF). We confirmed that in¬ 
creasing Ni) from 3 to 5 leaves the physics unchanged. 
This transition is independent of a smaller value of U2 
closer to that of ^^Rb, see Fig[ 2 ][c). For 0 > Oc, the con¬ 
densate density at the point (tt, 0 ) on the edge of the Bril- 
louin zone rises continuously as 0 increases further [Fig. [ 2 ] 
(a) and (b)]. Here, we observe the condensate density ris¬ 
ing continuously at (a^^(Q)), where Q = (tt, 0 ), for each 
rUz component [Fig. [ 2 ][a)], and acquires a non-zero imag¬ 
inary part in both momentum components (q = 0,Q), 
[inset of Fig. [ 2 ] (a)]. The precise location of this phase 
transition is determined by an interplay between the sin¬ 
gle particle spin-orbit physics and the underlying lattice, 
as we explain below. The non-zero momentum conden¬ 
sate is pinned to an ordering vector Qor commensurate 
with the underlying lattice Qor = Q for Oc < 0 < tt, 
as shown in Fig. [ 2 ][b). No other momentum components 
condense. 

In the continuum, the SDW phase results purely from 
non-interacting bosons condensing into multiple minima 
of the single-particle dispersio n^^i^^ . Here, the lattice 
suppresses condensation into incommensurate spin-orbit 
minima, until the spin-orbit wave vector is large enough 
such that the single-particle minima becomes close to 
(tt, 0 ). This is because in the limit of small on-site occu¬ 
pation, it costs a large amount of energy to condense 
away from a lattice site, as a result the superfluid is 


pinned to momenta commensurate with the lattice wave- 
vector rather than the spin-orbit wave-vector. We have 
checked that our calculations reproduce continuum re¬ 
sults at large on-site occupation (A'b ^ 10) and t/Uo^ 
where the lattice no longer plays a dominant role in the 
physics. 

A second non-trivial feature of the spin-1 system (com¬ 
pared to its pseudo spin-1/2 counterpart) is the appear¬ 
ance of a longitudinal SDW, concomitant with density 
wave order—. We take ( 5 '^(Q)) as the SDW order pa¬ 
rameter, which rises continuously upon entering the S-SF 
phase as shown in Fig. [2] (c). We find that the SDW is 
well described by {S^{x,y)) = As^ P Bs^ cos( 7 rx), where 
a = x^y^z and an amplitude Bs<^ on the order of the 
condensate density [see Fig. [ 2 ](d)]. 

To understand the origin of the longitudinal SDW, 
note that deep in the S-SF phase we can treat the spin 
order parameter at the mean field level, i.e. (S(q)) ^ 
Ek(«a(q + c,/3(a/3(k)). Therefore, the amplitude 

and wave-vector of the spin density wave order is com¬ 
pletely tied to the existence of a condensate at q = (0, 0) 
and (tt, 0 ). We also find a very weak charge density wave 
(not shown), however it is strongiy penaiized by the iat- 
tice and interactions. A major advantage of the spin-1 
system over the spin-1/2 counterpart is that the existence 
of additionai iongitudinal SDW order shouid make the 
stripe superfiuid phase easier to detect experimentaliy. 
In the spin-1/2 case with NIST type SOC, there is no 
iongitudinai SDW at zero detuning, and the contrast of 
the density osciilations is typicaiiy extremeiy smail^“— . 
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FIG. 3: Evolution of the S-SF as a function of the strength 
of interaction t/Fo for a fixed value of /x/Fo = 0.23, with 
translational invariance in the y-direction. Lx — U 2 — 
—0. 251 / 0 , = O.ll/o, ^ = TT, and Nh = 4. The SF order 

parameter for each nriz value at q = (0, 0) (a), and q = (tt, 0) 
(b), with labels for rUz shared between the two. The total 
spin (c) and the SDW order parameter (d) for both x (green) 
and z (cyan) components. The order parameters (*S'“(Q)) and 
(am;3(0,0)) are continuously suppressed at the S-SF to Q-SF 
transition [(a) and (d)], while the continuous Mott transition 
is captured by {S^) and (am^(Q)) going to zero [(b) and (c)]. 


Strong Coupling Phases — Having established the 
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weak coupling physics and the nature of the S-SF phase, 
we turn to the central question of this paper: how do 
strong interactions compete/complement the spin-orbit 
coupling to affect the condensate and spin/charge density 
wave orders? To address this, we now explore the phase 
diagram near the Mott-superfluid boundary at fixed SOC 
wave-vector = tt, and tune the ratios ji/Uo and t/Uo^ 
which will eventually drive the model to a Mott insulat¬ 
ing stated"—. 

Fig.[T]shows the resulting phase diagram at 0 = tt. We 
present results for Lx = 20 (with translational symmetry 
along y) with Nij = 4 and U 2 /UQ = —0.25. For large t/Uo 
we find the S-SF phase, which is well described by our 
previous discussion for a large region of the parameter 
space. As shown in Figs. [3] (a) and (d), decreasing t/Uo 
sends (a^^(q = 0)) —> 0, which destroys the SDW order 
continuously (5'"(Q)) ^ 0, giving way to a superfluid 
where (a^_^(Q)) 7^ 0 and ( 5 '"(q = 0)) 7^ 0 (a = x,z) 
[see Figs. [3] (b) and (c)], which we define as the Q-SF 
phase due to the condensate only appearing at a finite 
momentum. Increasing LIr first shifts the S-SF/Q-SF 
phase boundary to larger t, and finally completely polar¬ 
izes the gas. Decreasing 6>, shrinks the regime of stability 
of the Q-SF phase, nonetheless a 0 ^ tt is experimentally 
realistic. 

The phase transition from the S-SF to the Q-SF phase 
is a non-trivial interplay between competing energy scales 
in the problem. Upon decreasing t/Uo, interactions sup¬ 
press the overall amplitude of the condensate, and also 
favor condensation at a single momentum, because of the 
additional exchange energy between condensed bosons 
at q = (0,0) and (7r,0). In the continuum, local interac¬ 
tions induce a phase transition to a “plane wave state”—, 
where bosons either condense at one or the other single¬ 
particle minima, a spontaneously broken symmetry. In a 
lattice however, a state with a spatially modulating order 
parameter ((a^^(Q)) 7^ 0), but a homogeneous density 
has lower spin-orbit kinetic energy than a state with a 
uniform order parameter ((a^^(q = 0)) 7^ 0) and den¬ 
sity. This can be seen for ^ = tt by applying a spin 
rotation about the ^-axis {Sx —Sz and Sz Sx), 
which results in the hopping matrix element switching 
sign t ^ —t for ruz = ±1, favoring condensation at finite, 
rather than zero momentunt^. The Q-SF phase, where 
the superfluid order parameter has the same amplitude, 
but switches sign from site to site, minimizes the kinetic, 
interaction and Zeeman energies. This phase is therefore 
universally picked as the lower energy state. We remark 
that such finite momentum condensed phases of bosons 
are rather non-generic, and are actively being studied 
for bosons in ladder geometries with artificial magnetic 
fields^. 

As shown in Figs. [3][b) and (c), as t/Uo further de¬ 
creases {Sx)^ {am^ (Q)) vanishes continuously at the Mott 
transition, see Fig. [3l[b) and (d). Putting all of this to¬ 
gether we arrive at our main results shown in Fig. [TJ 
Importantly, the transition into each ferromagnetic Mott 
lobe is distinct from the case 0 = 0 (absence of SOC), 


since it is now the Q-SF phase which is becoming gapped 
out. Within the GMFT, each Mott lobe is homogeneous 
and ferromagnetic, as the Mott state is only captured at 
the t = 0 level and particle-hole and spin fluctuations are 
neglected. 



FIG. 4: Spin density wave of the S-SF phase in a harmonic 
trap with t = 0.2Uo, y = l.OUo, A = O.OlUo, Ui? = O.lUo, 
Lx — Ly — 35, 0 = 7T, and Nb = 5. (a) SDW for U 2 /U 0 = 
—0.25, the color displays the value of \{S^{x, y)) \ in the x — y 
plane, (b) Integrated spin density {S^{x)) = {x,y))/L 

as a function of x for U 2 /U 0 = —0.005 and —0.25. 

The excitation spectrum of the superfluid phases we 
find will be quite interesting. In the Q-SF, we expect a 
quadratically dispersing ferromagnetic spin mode about 
q = (0,0) and a density mode which is linearly dispers¬ 
ing at small |q| about Q, but has a roton-maxon like 
structure at intermediate q. This structure arises because 
bosons condense into only one of the single particle min¬ 
ima, effectively gapping out the other. The transition 
from the Q-SF to the S-SF should be marked by the van¬ 
ishing of the roton gap. A detailed calculation of the 
excitation spectrum in the spin-1 SOG Bose gas will be 
the subject of future work. 

Experimental Relevance — Our system can be read¬ 
ily studied in experiments by generalizing the existing 
setup of Gampbell et al. to an optical lattice^. Experi¬ 
mentally, the detection of {S^) is much more challenging 
than (5'^), while a precise measurement of each {umj is 
possible. Within the experimental apparatus of the NIST 
setup^, and the measurement basis we choose, the spin 
density wave in [Figs. [3] (c) and (d)] should be experi¬ 
mentally accessible. We have also numerically confirmed 
the existence of the S-SF phase for U2 = —O.OObUo ap¬ 
plicable to ^^Rb [see Figs. He) and [H[b)]. Our results 
in Figs. [T][3] provide a theoretical platform to interpret 
and understand experimental data. A key ingredient in 
most experiments is a harmonic trap Wi = X{x‘f -h y‘f) 
(we assume strong confinement in the 2)-direction, and 
use open boundary conditions). In Fig. [U^a) we plot the 
{S^{x^y)) order parameter for U2/U0 = —0.25 taking 
into account the external trapping profile, which shows 
oscillations that would be measurable in situ. Fig. [4l[b) 
shows the dependence of the magnitude of spin modula¬ 
tions within the trap on U 2 , where for a small value (ap¬ 
plicable to ^^Rb) the SDW has non-zero yet small mod¬ 
ulations, which increase dramatically for U2 = — 0.25Uo 
(closer to that of ^Li). Both the Q-SF and the S-SF can 
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also be readily studied in time-of-flight, which should re¬ 
veal a single peak for the former and a bimodal structure 
for the latter. Furthermore, the different excitation spec¬ 
tra in the two cases is another experimental probe, which 
can be measured using Bragg spectroscopy^. Our stud¬ 
ies highlight the non-trivial interplay between large spin, 
lattice physics and interactions and will serve to guide 
ongoing experiments on these systems. 
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